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ABSTRACT 

> ' 

• i-H , 

We have performed optical imaging observations of the dark cloud L1251 at 
multiple wavelengths, B, V, R, and I, using the 105 cm Schmidt telescope at 
the Kiso Observatory, Japan. The cloud has a cometary shape with a dense 
"head" showing star formation activity and a relatively diffuse "tail" without 
any signs of star formation. We derived extinction maps of As and Ay with a 
star count method, and also revealed the color excess (E B _ V , E V _ R , and E V -i) 
distributions. On the basis of the color excess measurements we derived the 
distribution of the ratio of total to selective extinction Ry over the cloud using 
an empirical relation between Ry and A\j Ay reported by Cardelli et al. In the 
tail of the cloud, Ry has a uniform value of ~3.2, close to that often found in 
the diffuse interstellar medium (~3.1), while higher values of Ry=A— 6 are found 
in the dense head. Since Ry is closely related to the size of dust grains, the high 
.Ry-values are most likely to represent the growth of dust grains in the dense 
star-forming head of the cloud. 



Subject headings: ISM: clouds — dust, extinction — ISM: individual (L1251) - 
stars: formation 
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1. Introduction 

The ratio of total to selective extinction, Ry=Ay / 'Eb-v ', is a measure of wavelength 
dependence of the interstellar extinction, representing the shape of the extinction curve . Ry 
is an important parameter to characterize dust properties at optical wavelengths because it is 
closely related to the size and composition of dust grains (e.g., Mathis, Rumpl, & Nordsieck 
1977; Hong & Greenberg 1978; Kim, Martin, & Hendry 1994). It is well established that 
R v has a uniform value of ~3.1 in the diffuse interstellar medium (e.g., Whittet 1992), 
but it tends to be much higher in dense molecular clouds (Ry=4— 6; e.g., Vrba, Coyne, & 
Tapia 1993; Strafella et al. 2001). Many theoretical predictions attribute the higher Ry 
in the dense cloud interior to a significant change of dust properties, particularly in size 
(e.g., Mathis 1990; Ossenkopf & Henning 1994). Since the timescale of grain growth due to 
coagulation of smaller grains and/or accretion of molecules in the gas phase onto the larger 
grain surface depends greatly on density (e.g., Tielens 1989), the higher Ry-values observed 
in dense molecular clouds are most likely to represent the growth of dust grains. 

There have been a number of observational reports on Ry measurements to date (e.g., 
Larson et al. 2000; Whittet et al. 2001). However, they are mostly based on the spec- 
trophotometry of a limited number of stars lying in the background of the clouds, and the 
-Ry-values have been measured only toward the selected stars. A large-scale distribution of 
Ry over an entire cloud, as well as its local variation inside, are therefore still poorly known. 

In order to investigate how much Ry can vary in a single cloud, we have carried out 
imaging observations of L1251 at multiple wavelengths, B, V, R, and /, using the 105 cm 
Schmidt telescope at the Kiso Observatory, Japan. This cloud is located in the Cepheus Flare 
(a~22 A 35"\ 5~75°00'; J2000.0) and has an elongated shape. Since L1251 is an isolated dark 
cloud located at a relatively high Galactic latitude (~15°), it is expected that the cloud is 
free from contamination due to unrelated clouds in the same line of sight. Kun & Prusti 
(1993) estimated its distance to be 300±50 pc with a color excess measurement and spectral 
classification of stars. They also found 12 H a emission-line stars at the eastern part of the 
cloud, five of which are detected as IRAS point sources. Global molecular distribution 
of the cloud was first revealed by Sato & Fukui (1989) and Sato et al. (1994) through 
13 CO and C 18 (J=l— 0) observations carried out with the 4 m radio telescope at the 
Nagoya University (HPBW=2.'7). Their observations clearly revealed the cometary "head- 
tail" morphology of L1251. They also found that two of the IRAS point sources located in 
the eastern "head" part of the cloud are associated with molecular outflows, while no sign of 
ongoing star formation was found in the western "tail." These characteristics of L1251 make 
this cloud a suitable site to investigate the possible difference in the dust properties (e.g., 
Ry) between star-forming and non-star-forming regions in a single cloud. 

In this paper, we report results of the BVRI observations toward L1251. The observa- 
tional procedures are described in § 2. The distribution of extinction (Ab and Ay) and color 
excess at each band (Eb~v, Ey^R, Ey_i) are derived from the obtained data, which reveal 
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the global dust distribution of the cloud. These are shown in §§ 3 and 4. We also derive 
Ry distribution over the cloud on the basis of color excess measurements and an empirical 
relation between Ry and A\j Ay (Cardelli, Clayton, & Mathis 1989). A significant difference 
has been found in Ry distribution between the dense head and the relatively diffuse tail of 
the cloud, suggesting a larger dust size in the dense head than in the tail. We introduce our 
procedure to create Ry map and compare the obtained Ry distribution with the other data 
in § 5. Our conclusions are summarized in § 6. 

2. Observations and Data Reduction 

We carried out imaging observations of L1251 from 2000 November 27 to December 2 
using the 105 cm Schmidt telescope equipped with the 2kCCD camera at the Kiso Observa- 
tory. The camera has a backside-illuminated 2048 x 2048 pixel chip covering a wide field of 
view of ~50'x50' with a scale of l."5 pixel -1 at the focal plane of the telescope. Using the 
camera, we mapped L1251 with the B, V, R, and / filters at a 40' grid spacing along the 
equatorial coordinates, so that adjacent frames could have an overlap of ~10' on the sky. 
We performed integrations of 200, 100, 100, and 50 s for the images obtained with the B, 
V, R, and / filters, respectively. The typical seeing was 3."0— 4."5 during the observations. 
In total, we obtained 10 frames per filter, covering the entire extent of the cloud. 

To calibrate the images, we obtained dome flat frames in the beginning of the observa- 
tions and frequently obtained bias frames during the observations. In addition, we observed 
some sets of standard stars listed by Landolt (1992) to correct the atmospheric extinction 
and calibration. 

We processed the calibrated images in a standard way using the IRAF packages. 1 After 
replacing values at the defected pixels in the detector (i.e., dead pixels) by those interpolated 
using surrounding pixels, we applied the stellar photometry package DAOPHOT (Stetson 
1987) to the images to detect stars and then to make photometric measurements. To im- 
prove the photometry we further measured the instrumental magnitudes of the detected 
stars by fitting a point-spread function (PSF) and then calibrated them for the atmospheric 
extinction. Among all of the detected stars we selected the ones having an intensity greater 
than 4.0 a noise level of the sky background. Some false detection, however, could not be 
completely ruled out, which we later removed by eye inspection. 

For calculating the plate solution, we compared the pixel coordinates of a number of 
the detected stars with the celestial coordinates of their counterparts in the Digitized Sky 



1 IRAF is distributed by the National Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under coorperative agreement with the National Science 
Foundation. 
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Survey (DSS), which resulted in a positional error of less than 1" (rms). All of the stars 
detected in our observations were summarized separately with respect to the filters we used 
(BVRI), which we will adopt to generate an extinction map using a star count method. 
The limiting magnitudes are 20.0, 19.7, 19.0, and 18.3 mag for the B, V, R, and / bands, 
respectively. Among all of the identified stars ~33,000 were detected in all of the four 
bands. These stars are useful for probing the distribution of the color excess in the cloud. 
We converted the instrumental magnitudes of these stars into the magnitudes in the standard 
Johnson-Cousins photometric system (Johnson & Morgan 1953; Cousins 1976), and we also 
measured their colors (B—V, V—R, and V—I). For the color excess measurements, we 
constructed a 5\^i?/-matched star list for ~30,600 stars whose photometric errors in the 
color measurements, <Jv-\, are less than 0.15 mag. One should note that we excluded all of 
the known young stellar objects (Sato & Fukui 1989; Sato et al.1994; Kun & Prusti 1993) 
cross-identified in our photometries in order to avoid errors in the following analyses since 
they are unrelated to the global extinction in L1251. 



3. Distribution of Extinction 

On the basis of the stars detected in the B and V band, we derived extinction (i.e., Ab 
and Ay) maps of L1251 using a star count method. Given a cumulative stellar density in 
the extinction- free reference field iV re f as a function of the A-band magnitude m x (i.e., the 
Wolf diagram; Wolf 1923), a logarithmic cumulative stellar density log N measured at the 
threshold magnitude m Xt0 can be converted into the extinction A\ as 

A x (l,b,m X fi) = m Xfi - f- 1 [\ogN(l,b,m x ,o)} , (1) 

where / and b are the Galactic coordinates, and / -1 is the inverse function of / defined as 
log iV re f = f(m x ). The function / is often assumed to be linear as / = a + bm x , where 
b= [dlogN rei /dm x } mx=mxo (e.g., Dickman 1978). 

To derive N(l,b,m Xt0 ) in equation (1) from the obtained data, we first set circular cells 
with a 4' diameter spaced by 2' along the Galactic coordinates, and we counted the number 
of stars in the cells, setting the threshold magnitude to be m B $ = 19.5 mag and m v $ = 
19.2 mag, which is well above the detection limit in our observations. We then smoothed 
the stellar density map to a 6' resolution with Gaussian filtering for a better signal-to-noise 
ratio to produce the final stellar density map N(l,b,m Xt o), which we adopted to derive the 
Ab and Ay distribution. After smoothing typical stellar density in the reference field and 
in the most opaque region was ~20 cell -1 and ~1 cell -1 for B band and ~30 cell -1 and ~2 
cell -1 for V band, respectively. 

We composed the stellar luminosity function f Te f(m x ) in the reference field, i.e., the area 
expected to be free from extinction in the observed region (see Fig. 1). We found that the 
function is actually not linear, but has a curved shape. We also found that stellar densities 
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at the threshold magnitude m A ,o systematically decrease along with the Galactic latitudes b. 
To deal with this problem, we fitted the average luminosity function f Te {(m\) by a 4th-order 
polynomial and scaled it to the background stellar density iV back at m A =m A as 

/(i,^^ l08 y, M - m>j) t^;, (2) 

Jref {m\,0) ^ 

where c n represents the best-fitting coefficients. The value f r cf( m x,o) corresponds to the mean 
stellar density at m\$ in the reference field. We determined iVback^, b, m\ t o) by assuming an 
exponential dependence on b as 

AW (l,b,m Xt o) = iVoexp (-a \b\) , (3) 

where Nq and a are constants determined using values of N(l, b, mj^) in the reference field. 
We obtained the coefficients (N , a) to be (61.5, 0.07) for B band and (106.2, 0.09) for V 
band. 

We numerically calculated the inverse function of / in equation (2), and used it in 
equation (1) to derive As and Ay after substituting the actually measured stellar density 
N(l,b,mxfi) into The resulting Ay map is shown in Figure 1. The distribution of As 
(not shown) is very similar to that of Ay. Total uncertainty in A\ varies from region to 
region, depending on the number of stars falling in the cells. Typical uncertainties in the 
reference field (A\ ~ mag) and in the most opaque region are estimated to be ~0.5 and 
~0.8 mag for the B band, and ~0.35 and ~0.7 mag for the V band, respectively. 

The Ay map reveals an interesting dust distribution in L1251 with a cometary head-tail 
morphology, which is similar to the molecular distribution traced in 13 CO (Sato et al. 1994). 

4. Distribution of Color Excess 

We derived the color excess distribution in L1251 using the near-infrared color excess 
(NICE) method originally developed by Lada et al. (1994). The method compares the 
color index of stars in the reference field with those in the dark clouds. If we assume that 
the population of the stars across the observed field is invariable, the mean stellar color 
in the reference field can be used to approximate the mean intrinsic color of stars (i.e., 
<V— \> int =<V— A> rcf ), and the color excess can be derived by subtracting <V—\> Ici from 
the observed color of stars. To obtain the color excess maps, we adopted an "adaptive grid" 
technique (Cambresy 1998, 1999) to the ~30,600 stars detected in all of the four bands (B, 
V, R, and /). This technique consists in fixing the number of counted stars N per variable 
cell to achieve a uniform noise level in the resulting map. To estimate the mean color excess 
<Ey^\> (X=B, R, or /) in circular cells spaced by 2' along the Galactic coordinates with 
undefined sizes, we calculated the mean color of the closest N stars to the center positions of 



- 6- 



the cells by setting iV to be 10 and derived their shifts from the mean color in the reference 
field <V— A> rcf as 

(l,b)-(V- A) ref , (4) 

where (V— A)j is the color index of the ith star in a cell. The distance to the Nth star from 
the cell center gives the radius of the cell, corresponding to one half of the angular resolution. 
<V— A> rc f is the mean color of the ~2000 cells in the reference field, and we determined 
<B—V> TC f ) <V— R> ie ii and <V— 7> re f to be 1.02, 0.55, and 1.23 mag, respectively. 

In the reference field, we found that <V— A> slightly but systematically decreases with 
increasing Galactic latitudes b. This may be caused by the presence of diffuse dust layer of 
the Galaxy and/or gradual change in the stellar population toward the Galactic plane. We 
calibrated the effect by substituting <V— A> ref in equation (4) for the following equation 
assuming an exponential dependence on b as 

(V- A) (l,b) = Aexp(-a|6|) , (5) 

where A and a are the fitting coefficients. We determined the coefficients (A, a) to be (1.60, 
0.03), (2.60, 0.05), and (0.84, 0.03) for <B-V> rcf , <V-R> rcf , and <V-I> rci , respectively. 

We then smoothed the color excess maps with a Gaussian filter (FWHM=4') to reduce 
noise. The resulting resolution of the maps is typically ~5' in the reference field and ~10' in 
the most opaque region. An example of the obtained color excess maps is shown in Figure 2, 
where we show the spatial distribution of Eb-v- The distributions of Ey-n and Ey-i (not 
shown) are similar to that of Eb-v- L1251 is also traced in the all-sky Eb-v map with 6' 
resolution (Schlegel, Finkbeiner, & Davis 1998), which is derived from far- infrared emission 
data (COBE/DIRBE, IRAS/ISSA). The cloud appears similar in both the all-sky map 
and our color excess maps, indicating that the optical color excess and the far-infrared 
thermal emission of dust trace the same region. 

Uncertainty in deriving <Ey-\> is given by following equation (e.g., Arce & Goodman 
1999): 

1/2 

(6) 

where a(y-\) i is the photometric error in V— A of the ith star in each cell. The value cr ce ii 
is the standard deviation of the iV-sampling distribution of the mean, where the original 
distribution is the V— A color of the stars in the reference field. It is equivalent to the 
standard deviation of the <Ey-\> distribution in the reference field. We estimated a cc n to 
be 0.04, 0.03, and 0.05 mag for the color of B—V, V—R, and V—I, respectively. In the 
E B -v, E V -r, and E V -i maps, a typical a <Ev _ x> -va\ue calculated using equation (6) is 0.05, 
0.03, and 0.05 mag, respectively. 
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We derived the color excess diagram E B -v versus E V -i to see the global trend of 
reddening by dusts in L1251, which is shown in Figure 3. We investigated the relation 
separately in the star-forming head (7>114?2) and in the non-star-forming tail (/<114?2). 
For comparison, we also show Eb-v versus Ey_j relations in the cases of Ry =3.1 and 6.5, 
adopting an empirical extinction law reported by Cardelli et al. (1989). 

It is clear that there is a significant difference in the Eb-v versus Ey_j relations obtained 
in the two parts of the cloud. In the non-star-forming tail the data are distributed along 
the line of Ry—3.1, a typical value in the diffuse ISM, and the linear least-squares fit to the 
data actually yields a slope of Ey-i/ ^B-y=1.21, which corresponds to Ry=3.2. On the 
other hand, the data in the star-forming head are distributed above the line of Ry =3.1 in 
the range of E B -y>0A mag, and they gradually get close to the line of Ry=6.5 for larger 
E B -v- It is evident that optical properties of dust grains significantly vary even in a single 
cloud. The different extinction laws are well separated into two groups only when we divide 
them into star-forming and non-star-forming parts. Since Ry increases with increasing size 
of dust grains, larger Ry- value in the star- forming head is indicative of the evolution (i.e., 
growth in size) of dust grains. We note that Ey-i increases in the same way (along with 
i?y~3.1) up to _Eb-v~0.4 mag both in the tail and in the head part, which indicates dust 
optical properties remain normal at the diffuse boundary of the cloud. 



5. Derivation of R v Distribution and Evidence for Grain Growth 

In order to investigate the local variation of the dust properties in L1251, we derived 
the Ry distribution using a relationship between Ry and Ey_\j Eb-v that we produced 
from the i?y-dependent extinction law derived by Cardelli et al. (1989). They found that 
observed A\j Ay versus R v l relations toward a number of line-of-sight stars are well fitted 
over the wavelength range from 0.125 /im to 3.5 fim as 

(A x /A v ) = a(x) + b^Ry 1 , (7) 

where x is A -1 jum -1 . The wavelength-dependent coefficients a(x) and b(x) are shown in 
equations (2a-2b) and (3a-3b) of Cardelli et al. (1989). By using equation (7) we reproduced 
Ry and E V -x/ 'Eb-v relations with the color excess ratios as E V -x/ E B -v=(A v -A\) /(A B - 
Ay). 

Color excess ratios of E V _ R / E B -v an d E v _i/E B -v are used to obtain the i?y-values at 
the cells that we used to derive the color excess maps. In each cell we numerically calculated 
the iiV-value, which minimizes the sum of residual variance; i.e., \{Ev-r/ Eb-v) ob s — {Ev-r/ 
E B -v)xnod c \} 2 + {{E V -i I ' E B -v)obs-{E V -i I 'EB-v)modc\] 2 '■ We carried out the Ry fitting only 
in the region with i?^_y>0.4 mag in order to avoid a large error in Ry arising from too 
small Eb-v values. The resulting Ry map is shown in Figure 4. As seen in the figure, Ry 
in the star-forming head is systematically higher than in the tail with the maximum value 
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of i?y~6. The mean value of Ry in the head and tail are ~4.6 and ~3.2, respectively. The 
feature in the Ry distribution is consistent with what we found in the color excess diagrams 
derived in the previous section (Fig. 3). 

Uncertainties in the Ry measurement, (Jr v , are estimated using Monte Carlo simu- 
lations. We first generated noise that follows a Gaussian distribution with the standard 
deviation equal to the 1 a error of Ey_\/ Eb-v measurements, and then we performed the 
Ry fitting to the noise-added Ey_\/ Eb-v data. We repeated this calculation 1000 times 
in each cell, and regarded the standard deviation of the resulting i?y-values as <Jr v . The 
value of <jr v varies from ~0.5 to ~0.9 (0.69 on the average) and tends to be relatively higher 
in cells with small Eb~v- In addition to the random noise our Ry map may suffer from 
a systematic error due to a number of unidentified pre-main-sequence (PMS) stars formed 
in the dense head. A certain fraction of PMS stars are known to have a blue color excess 
(e.g., Kenyon & Hartmann 1995), which may result in an overestimate in Ry. However, this 
possible error should be negligible in our case, as we discuss in the Appendix. 

We further checked the validity of the color-excess-based R v measurements by using A B 
and Ay derived with star count method. Though large uncertainties (typically ~0.5 mag in 
our data) are included in Ab and Ay data, Ry can be directly computed as Ry=Ay / (Ab — 
Ay) without using an empirical model. We made Ay versus Ab diagram for the same 
cells used in Ry measurements and separated them into the star- forming head (/>114?2) 
and non-star-forming tail (/<114?2). The slope of the Ay versus Ab relations yields mean 
i?y-value of data points as <Ry>=l/(a— 1), where a=AB/Ay. The linear least-squares 
fit to each dataset resulted in ah ea d=l-23±0.09 and a ta ii=l-34±0.07, which correspond to 
<i?y>head~4.3 and <i?v r >taii~3.0, respectively. These values are consistent with the results 
of color-excess-based Ry measurements. 

In Figure 4 it is particularly noteworthy that there is a good coincidence between the 
local peaks of the Ry distribution and the locations of YSOs. There are three strong Ry 
peaks in the head at (/,6) = (114?60, 14?47; R v =6.1), (114?50, 14?67; R v =5.8), and (114?27, 
14?77; R V =5.Q). The first two peaks show a good agreement with the locations of the two 
outflow sources reported by Sato & Fukui (1989) and Sato et al. (1994), though the third 
peak has no apparent corresponding source. In addition, the three Ry peaks also coincide 
with the positions of molecular cloud cores found in C 18 (named "B", "C", and "E" by Sato 
et al. 1994; see their Table 1 and Fig. 2), indicating that the region with high Ry- values is 
rich in molecular gas. 

To investigate the dependence of Ry on the dust column density, we derived a diagram 
of Ay versus Ry. For a direct comparison we remeasured the Ay-values with the same cells 
used to produce the Ry map. In Figure 5 we show the diagram measured both in the dense 
head and in the diffuse tail. It is clear that the Ry increases along with Ay in the head, 
while that in the tail remains constant regardless of Ay. 
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All of the above results obtained from the comparison of the Ry map with the other data, 
i.e., the coincidence of the Ry peaks and the outflow sources, as well as the C 18 cores, and 
the difference in the Ay versus Ry relation between the head and the tail, strongly indicate 
that Ry is higher in the dense star-forming head than those in the diffuse tail without any 
signs of star formation. Because a higher value of Ry is expected for dusts with a larger size, 
the dusts in the head of L1251 may actually be larger than in the diffuse tail, indicating the 
growth of dust grains through coagulation and accretion of smaller dust particles as well as 
molecules in gas phase onto larger dust grains in the dense molecular cloud cores experiencing 
star formation. The Ry distribution derived in this study may not be sufficient to probe 
the densest region of L1251 (A V >A— 5 mag) in detail because of a large optical depth at the 
observed wavelengths. However, a similar change of dust properties most likely representing a 
large dust size in molecular clouds was recently reported through submillimeter observations 
of the dust emission (e.g., MCLD 123.5+24.9: Bernard et al. 1999; Taurus: Stepnik et al. 
2003). The idea of grain growth is also supported by their study based on the optically much 
thinner submillimeter observation, implying that dust grains with a larger size than in the 
diffuse ISM may be common in dense molecular clouds. 

Our procedure to derive the Ry distribution on the basis of the average color excess 
measurements is useful to investigate a large-scale variation of the dust optical properties. 
It may be applicable also to the near infrared data in the future, with which we should be 
able to measure the Ry distribution in the densest region of L1251 to confirm the significant 
change of Ry therein. 

6. Summary 

We have carried out optical observations toward L1251 at multiple bands, B, V, R, 
and /, using the 105 cm Schmidt telescope equipped with the 2kCCD camera at the Kiso 
Observatory in Japan. Extinction and color excess maps of the cloud were derived using a 
star count method and a NICE method. From the observed E v _ x / E B -v (A=-R,-0, we also 
derived the Ry distribution using an empirical Ry versus A\j Ay relation (Cardelli et al. 
1989). 

(1) We investigated the color excess distribution of the cloud on the basis of the ~30,600 
stars detected in our observations. Color excess diagram, E B _ V versus E v _j, indicates a 
significant difference in Ry between the star-forming dense head (7>114?2) and the non- 
star-forming diffuse tail (K114?2). 

(2) We derived the Ry distribution over the cloud on the basis of the color excess 
measurements and an empirical relation between Ry and A x j Ay (Cardelli et al. 1989). From 
the resulting Ry map we found that Ry is systematically higher in the dense head (Ry=4— 6) 
than in the diffuse tail (i?y~3.2), which is consistent with what we found in the E B -v 
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versus E V -i diagram. We confirmed the validity of color-excess-based Ry measurements by 
comparing them with independently measured iiV-values; i.e., A V /(A B —A V ), which were 
directly computed using A B and Ay data from star count analysis. 

Local peaks of the Ry distribution in the head of L1251 coincide with the positions of 
the young stars accompanied by molecular outflows and the molecular cloud cores identified 
in C 18 0. We also found that Ry increases with increasing Ay in the dense star-forming head, 
while Ry is likely to be independent of Ay in the diffuse tail. Since Ry is closely related to 
the size of dusts, we suggest that these results are indicative of the growth of dust grains 
in the dense head due to coagulation and accretion of smaller dust particles and molecules 
onto the larger dust grains. 

We are very grateful to M. Tamura, K. Tatematsu, and T. Mizuno for their helpful 
comments and suggestions. Thanks are due to the staff at Kiso Observatory including K. 
Nakada, T. Miyata, T. Soyano, K. Tarusawa, Y. Tanaka, H. Mito, T. Aoki, N. Ito, and S. 
Nishiura for their kind support during the observations. This work was financially supported 
by the Grant-in- Aid for Scientific Research by the Japanese Ministry of Education, Science, 
Sports and Culture (Nos. 13640233 and 14022214). 

A. A Possible Systematic Uncertainty in Ry 

Here we consider the degree of uncertainties caused by the presence of unknown faint 
PMS stars in the star-forming head of L1251, which escaped detection in the previous H a 
survey (Kun & Prusti 1993). It is well known that a certain fraction of PMS stars show 
the optical blue excess, which is prominent in U—B (e.g., Kenyon & Hartmann 1995) and 
possibly in B—V. Their blue intrinsic color should result in an underestimate of Eb-v, for 
which we may overestimate Ry. 

In order to examine the fraction of blue excess PMS stars at B—V, we made a B—V 
versus V —I diagram of 96 PMS stars in the Taurus-Auriga region (see Table Al of Kenyon 
& Hartmann 1995). We found 23 (~25%) sources show a blue excess of greater than 0.2 mag 
compared with the main-sequence stars, while the rest of the sources have a similar color to 
the main-sequence stars with some reddening. 

If we assume that the ratio of the blue excess PMS stars measured in the Taurus- Auriga 
region is also valid in L1251, only two to three stars (25% of iV=10) in each counting cell are 
expected to be blue excess PMS stars even in an extreme case that all of the stars detected 
in our observations are PMS stars. However, at the maximum estimate of the blue excess 
PMS stars (about two to three stars per cell), there might be an influence on the value of 
E B ~v calculated as the average color of 10 stars (eq. [4]). We therefore recalculated E B -v 
taking the median color of 10 stars instead of the average, because it should be much less 
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affected by a few blue excess PMS stars. We compared the two values of E B -v, the median 
and the average colors of 10 stars, and found that they are quite consistent with each other 
even in the densest head of L1251 with a typical difference of only ~0.02 mag. This means 
that the possible error in Ry due to the blue excess PMS stars is negligible, and the high 
.Ry-value derived in the star-forming dense head is most likely to be real. 
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Galactic Longitude (degree) 

Fig. 1. — Ay distribution in L1251. The angular resolution is 6'. Contours start from 1 
mag with a step of 0.5 mag. Filled circles with white outlines represent the positions of the 
IRAS sources associated with molecular outflows (Sato et al. 1994). Plus signs represent 
the positions of the H a emission-line stars found by Kun & Prusti (1993). The observed area 
is delineated by the broken line. The area enclosed by the gray broken line is the reference 
field expected to be free from extinction. Typical uncertainties in Ay are ~0.35 mag and 
~0.7 mag, which are found in the reference field (Ay~0 mag) and in the most opaque region 
(Ay~5 mag), respectively. 
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Fig. 2. — Eb-v distribution in L1251. Contours start from 0.2 mag with a step of 0.1 mag. 
The uncertainty in E B _ V is ~0.05 mag. The highest and lowest angular resolutions achieved 
in the map are ~5' in the reference field and ~10' in the most opaque region, respectively. 
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Fig. 3. — (a) E B -v vs. E V -i relation for the non-star-forming tail of L1251 (/<114?2). The 
solid line denotes the linear least-squares fit (£'v-//-E , B-v r= l-21 corresponding to i?y=3.2) 
to all of the data points in the tail. The dotted and broken lines are for the cases of i?y=3.1 
and 6.5, respectively. The typical error in Eb-v and Ey^i is ~0.05 mag. (b) Same as (a) 
but for the star-forming head of L1251 (/>114?2). 
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Fig. 4. — R v distribution in L1251. The small dots represent the positions where we mea- 
sured the i?y-values. Filled circles with white outlines denote the IRAS sources associated 
with molecular outflows (Sato & Fukui 1989; Sato et al. 1994). Contours start from i?y=4 
with a step of 0.3. 
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Fig. 5. — (a) Ay vs. it!y relation for the non-star-forming tail of L1251 (7<114?2). The 
broken line denotes the mean Ry-value 3.2. A typical error in Ay and Ry is 0.5 mag and 
0.69, respectively, (b) Same as (a) but for the star-forming head of L1251 (/>114?2). 



